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Abstract. We establish a direct link between the time-reversal technique and 
the so-called adjoint method for imaging. Using this relationship, we derive 
new solution strategies for an inverse problem which arises in telecommunication. 
These strategies are all based on iterative time-reversal experiments, which try 
to solve the inverse problem experimentally instead of computationally. We will 
focus in particular on examples from underwater acoustic communication and 
wireless communication in a Multiple-Input Multiple-Output (MIMO) setup. 
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1. Introduction 

Time-reversal techniques have attracted great attention recently due to the large 
variety of interesting potential applications. The basic idea of time-reversal (often also 
referred to as 'phase-conjugation' in the frequency domain) can be roughly described 
as follows: A localized source emits a short pulse of acoustic, electromagnetic or 
elastic energy which propagates through a richly scattering environment. A receiver 
array records the arriving waves, typically as a long and complicated signal due to 
the complexity of the environment, and stores the time-series of measured signals 
in memory. After a short while, the receiver array sends the recorded signals time- 
reversed (i.e. first in - last out) back into the same medium. Due to the time- 
reversibility of the wave fields, the emitted energy backpropagates through the same 
environment, practically retracing the paths of the original signal, and refocuses, 
again as a short pulse, on the location where the source emitted the original pulse. 
This, certainly, is a slightly oversimplified description of the rather complex physics 
which is involved in the real time-reversal experiment. In practice, the quality of the 
refocused pulse depends on many parameters, as for example the randomness of the 
background medium, the size and location of the receiver array, temporal fluctuations 
of the environment, etc. Surprisingly, the quality of the refocused signal increases 
with increasing complexity of the background environment. This was observed 
experimentally by M. Fink and his group at the Laboratoire Ondes et Acoustique 
at Universite Paris VII in Paris by a series of laboratory experiments (see for example 
| 17l 1201 \21\). and by W. A. Kuperman and co-workers at the Scripps Institution 
of Oceanography at University of California, San Diego, by a series of experiments 
performed between 1996 and 2000 in a shallow ocean environment (see for example 

inillHlEEl) 

The list of possible applications of this time-reversal technique is long. In an 
iterated fashion, resembling the power method for finding maximal eigenvalues of a 
square matrix, the time-reversal method can be applied for focusing energy created 
by an ultrasound transducer array on strongly scattering objects in the region of 
interest. This can be used for example in lithotripsy for localizing and destroying 
gall-stones in an automatic way, or more generally in the application of medical 
imaging problems. Detection of cracks in the aeronautic industry, or of submarines 
in the ocean are other examples. See for example |39[ I40j . and for related work 
|7I I12|. The general idea of time-reversibility, and its use in imaging and detection, is 
certainly not that new. Looking into the literature for example of seismic imaging, the 
application of this basic idea can be found in a classical and very successful imaging 
strategy for detecting scattering interfaces in the earth, the so-called 'migration' 
technique However, the systematic use and investigation of the time-reversal 

phenomenon and its experimental realizations started more recently, and has been 
carried out during the last 10-15 years or so by different research groups. See 
for example [131111201121112312111121113 for experimental demonstrations, and 
Dlll3l3E10IiniII3l23l2HlEniEHll31lllE3llSlll3 for theoretical and 
numerical approaches. 

One very young and promising application of time-reversal is communication. 
In this paper we will mainly concentrate on that application, although the general 
results should carry over also to other applications as mentioned above. 

The paper is organized as follows. In section |21 we give a very short introduction 
into time-reversal in the ocean, with a special emphasis on underwater sound 
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communication. Wireless communication in a MIMO setup, our second main 
application in this paper, is briefly presented in section |3| In section 0] we discuss 
symmetric hyperbolic systems in the form needed here, and examples of such systems 
are given in section |S1 In section |S1 the basic spaces and operators necessary for 
our mathematical treatment are introduced. The inverse problem in communication, 
which we are focusing on in this paper, is then defined in section In section |H1 
we derive the basic iterative scheme for solving this inverse problem. Section El gives 
practical expressions for calculating the adjoint communication operator, which plays 
a key role in the iterative time-reversal schemes presented in this paper. In section 
1101 the acoustic time-reversal mirror is defined, which will provide the link between 
the 'acoustic time-reversal experiment' and the adjoint communication operator. The 
analogous results for the electromagnetic time-reversal mirror are discussed in section 
111! Section El combines the results of these two sections, and explicitly provides 
the link between time-reversal and the adjoint imaging method. Sections US El 
and 1151 propose then several different iterative time-reversal schemes for solving the 
inverse problem of communication, using this key relationship between time-reversal 
and the adjoint communication operator. The practically important issue of partial 
measurements (and generalized measurements) is treated in section lTHl Finally, section 
llTl summarizes the results of this paper, and points out some interesting future research 
directions. 

2. Time-reversal and communication in the ocean 

The ocean is a complex wave-guide for sound ^U] ■ In addition to scattering effects at 
the top and the bottom of the ocean, also its temperature profile and the corresponding 
refractive effects contribute to this wave-guiding property and allow acoustic energy to 
travel large distances. Typically, this propagation has a very complicated behaviour. 
For example, multipathing occurs if source and receiver of sound waves are far away 
from each other, since due to scattering and refraction there are many possible paths on 
which acoustic energy can travel between them. Surface waves and air bubbles at the 
top of the ocean, sound propagation through the rocks and sedimentary layers at the 
bottom of the ocean, and other effects further contribute to the complexity of sound 
propagation in the ocean. When a source (e.g. the base station of a communication 
system in the ocean) emits a short signal at a given location, the receiver (e.g. a user 
of this communication system) some distance away from the source typically receives 
a long and complicated signal due to the various influences of the ocean to this signal 
along the different connecting paths. If the base station wants to communicate with 
the user by sending a series of short signals, this complex response of the ocean to the 
signal needs to be resolved and taken into account. 

In a classical communication system, the base station which wants to 
communicate with a user broadcasts a series of short signals (e.g. a series of 'zeros' and 
'ones') into the environment. The hope is that the user will receive this message as a 
similar series of temporally well-resolved short signals which can be easily identified 
and decoded. However, this almost never occurs in a complex environment, due to 
the multipathing and the resulting delay-spread of the emitted signals. Typically, 
when the base station broadcasts a series of short signals into such an environment, 
intersymbol interference occurs at the user position due to the temporal overlap of the 
multipath contributions of these signals. In order to recover the individual signals, 
a significant amount of signal processing is necessary at the user side, and, most 
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importantly, the user needs to have some knowledge of the propagation behaviour 
of the signals in this environment (i.e. he needs to know the 'channel'). Intersymbol 
interference can in principle be avoided by adding a sufhcient delay between individual 
signals emitted at the base station which takes into account the delay-spread in the 
medium. That, however, slows down the communication, and reduces the capacity 
of the environment as a communication system. An additional drawback of simply 
broadcasting communication signals from the base station is the obvious lack of 
interception security. A different user who also knows the propagation behaviour 
of signals in the environment, can equally well resolve the series of signals arriving at 
his location and decode them. 

Several approaches have been suggested to circumvent the above mentioned 
drawbacks of communication in multiple-scattering environments. Some very 
promising techniques are based on the time-reversibility property of propagating wave- 
fields HSHiniliniEHMlESlEZllSl. The basic idea is as follows. The user who wants 
to communicate with the base station, starts the communication process by sending 
a short pilot signal through the environment. The base station receives this signal as 
a long and complex signal due to multipathing. It time-reverses the received signal 
and sends it back into the environment. The backpropagating waves will produce a 
complicated wave-field everywhere due to the many interfering parts of the emitted 
signal. However, due to the time-reversibility of the wave-fields, we expect that the 
interference will be constructive at the position of the user who sent the original pilot 
signal, and mainly destructive at all other positions. Therefore, the user will receive at 
his location a short signal very similar to (ideally, a time-reversed replica of) the pilot 
signal which he sent for starting the communication. All other users who might be in 
the environment at the same time will only receive noise speckle due to incoherently 
interfering contributions of the backpropagating field. If the base station sends the 
individual elements ('ones' and 'zeros') of the intended message in a phase-encoded 
form as a long overlapping string of signals, the superposition principle will ensure 
that, at the user position, this string of signals will appear as a series of short well- 
separated signals, each resembling some phase-shifted (and time-reversed) form of the 
pilot signal. 

In order to find out whether this theoretically predicted scenario actually takes 
place in a real multiple-scattering environment like the ocean, Kuperman et al have 
performed a series of four experiments in a shallow ocean environment between 1996 
and 2000, essentially following the above described scenario. The experiments have 
been performed at a Mediterranean location close to the Italian coast. (A similar setup 
was also used in an experiment performed in 2001 off the coast of New England which 
has been reported in Yang jlH]-) A schematic view of these experiments is shown 
in figure n The single 'user' is replaced here by a 'probe source', and the 'source- 
receive array' (SRA) plays the role of the 'base station'. An additional vertical receive 
array (VRA) was deployed at the position of the probe source in order to measure 
the temporal and spatial spread of the backpropagating fields in the neighbourhood of 
the probe source location. In this shallow environment (depths of about 100 — 200 m, 
and distances between 10 and 30 km) the multipathing of the waves is mostly caused 
by multiple reflections at the surface and the bottom of the ocean. The results of the 
experiments have been reported in jl8[ 1261 1ST] . They show that in fact a strong spatial 
and temporal focusing of the backpropagating waves occurs at the source position. A 
theoretical explanation of the temporal and spatial refocusing of time-reversed waves 
in random environments has been given in Blomgren et al 
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3. The MIMO setup in wireless communication 

The underwater acoustics scenario described above directly carries over to situations 
which might be more familiar to most of us, namely to the more and more 
popular wireless communication networks using mainly electromagnetic waves in the 
microwave regime. Starting from the everyday use of cell-phones, ranging to small 
wireless-operating local area networks (LAN) for computer systems or for private 
enterprise communication systems, exactly the same problems arise as in underwater 
communication. The typically employed microwaves of a wavelength at about 10- 
30 cm are heavily scattered by environmental objects like cars, fences, trees, doors, 
furniture etc. This causes a very complicated multipath structure of the signals 
received by users of such a communication system. Since bandwidths are limited and 
increasingly expensive, a need for more and more efficient communication systems is 
imminent. Recently, the idea of a so-called multiple-input multiple-output (MIMO) 
communication system has been introduced with the potential to increase the capacity 
and efficiency of wireless communication systems |22j. The idea is to replace a single 
antenna at the base station which is responsible for multiple users, or even a system 
where one user communicates only with his own dedicated base antenna, by a more 
general system where an array of multiple antennas at the base station is interacting 
simultaneously and in a complex way with multiple users. A schematic description of 
such a MIMO system (with seven base antennas and seven users) is given in figure [3 
See for example |^ 02 for recent overviews on MIMO technology. 

Time-reversal techniques are likely to play also here a key role in improving 
communication procedures and for optimizing the use of the limited resources 
(especially bandwidth) which are available for this technology 22.. 32., 31, ■ One big 
advantage of time-reversal techniques is that they are automatically adapted to the 
complex environment and that they can be very fast since they do not require heavy 
signal processing at the receiver or user side. In |^ , an iterated time-reversal 
scheme for the optimal refocusing of signals in such a MIMO communication system 
was proposed, which we will describe in more details in section IT^ 



Time-reversal and the adjoint method with an application in telecommunication 6 



users 



pilot signal a (t) 



o 



complex 
multiple 



base antennas 



scattering 



environment 



A[ o 

o 
o 

A4 o 

o 
o 



received signals 




Figure 2. Schematic view of a typical MIMO setup in wireless communication. 



In the present paper we will establish a direct link between the time- 
reversal technique and solution strategies for inverse problems. As an application 
of this relationship, we will derive iterative time-reversal schemes for the 
optimization of wireless or underwater acoustic MIMO communication systems. The 
derivation is performed completely in time-domain, for very general first order 
symmetric hyperbolic systems describing wave propagation phenomena in a complex 
environment. One of the schemes which we derive, in a certain sense the 'basic one', 
will turn out to be practically equivalent to the scheme introduced in |34[I35| . although 
the derivation uses different tools. Therefore, it provides a new interpretation of that 
scheme. The other schemes which we introduce are new in this application, and can be 
considered as either generalizations of the basic scheme, or as independent alternatives 
to that scheme. Each of them addresses slightly different objectives and has its own 
very specific characteristics. 

4. Symmetric hyperbolic systems 

We treat wave propagation in communication systems in the general framework of 
symmetric hyperbolic systems of the form 

u(x,0) = 0. (2) 

Here, u(x, <) and q(x, i) are real-valued time-dependent iV-vectors, x £ R^, and 
t £ [0,T]. r(x) is a real, symmetric, uniformly positive definite N x A'^-matrix, i.e., 
r(x) > e for some e > 0. Moreover, <I>(x) is a symmetric positive semi-definite N x N 
matrix, i.e., <&(x) > 0. It models possible energy loss through dissipation in the 
medium. The are real, symmetric and independent of (x, i). We will also use the 
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short notation 

i—l 

In addition to the above mentioned assumptions on the coefficients r(x) and $(x), we 
will assume throughout this paper that all quantities r(x), <&(x), u(x, i) and q(x, t) 
are 'sufficiently regular' in order to safely apply for example integration by parts and 
Green's formulas. For details see for example |ll[l33j . Throughout this paper we will 
assume that no energy reaches the boundaries dfl during the time [0, T], such that we 
will always have 

u(x,i)=0 on dnx[0,T]. (4) 
The energy density f (x, t) is defined by 

- i(r(x)u(x,t),u(x,i))^ 

1 ^ 

= 2 X! r„„(x)u„(x,<)u„(x,t). 

m^n—l 

The total energy £{t) in at a given time t is therefore 

£{t) = i^(r(x)u(x,t),u(x,i))^ dx. 
The Hux J-{x, t) is given by 

^,(x,i) = i(i?*u(x,t),u(x,i))^ , z = 1,2,3. 



5. Examples for symmetric hyperbolic systems 

In the following, we want to give some examples for symmetric hyperbolic systems as 
defined above. Of special interest for communication are the system of the linearized 
acoustic equations and the system of Maxwell's equations. We will discuss these two 
examples in detail in this paper. Another important example for symmetric hyperbolic 
systems is the system of elastic waves equations, which we will however leave out in 
our discussion for the sake of shortness. We only mention here that all the results 
derived here apply without restrictions also to linear elastic waves. Elastic wave 
propagation becomes important for example in ocean acoustic communication models 
which incorporate wave propagation through the sedimentary and rock layers at the 
bottom of the ocean. 



5.1. Linearized acoustic equations 

As a model for underwater sound propagation, we consider the following linearized 
form of the acoustic equations in an isotropic medium 

/o(x)-^ + gradp = q„ (5) 

+ divv = qp (6) 
p(x,0) = 0, v(x,0)=0. (7) 
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Here, v is the velocity, p the pressure, p the density, and k the compressibiUty. We 
have TV = 4, u = (v,p)^ and q — (q„, q^)^ (where 'T' as a superscript always means 
'transpose'). Moreover, wc have 

r(x) = diSL^{p{x) , p{x) , p{x) , k{x)) and $(x) = 0. 

With the notation ip = {81,82,83)^ , we can write A as 

The operators i = 1, 2, 3, can be recovered here from A by putting 

j^i ^ f 1 where Am,n = 9i, 
1^ elsewhere 

The energy density £{x,t) is given by 

£{x,t) = i(p(x)|v(x,i)p + K(x)/(x,t)), 

and the energy flux J^{x., t) is 

T{x,t) = p(x,t)v(x,t). 

We mention that the dissipativc case $(x) 7^ can be treated as well in our framework, 
and yields analogous results to those presented here. 



5.2. Maxwell's equations 

As a second example, wc will consider Maxwell's equations for an anisotropic medium 
with some energy loss due to the inherent conductivity. This can for example model 
wireless communication in a complex environment. 
8E 

e(x)— - curlH + a(x)E = q^ (8) 

m(x)— +curlE =ciH (9) 

E(x,0) = 0, H(x,0) = 0. (10) 

We have N = 6, and u = (E, H)^, q = (q^;, q^)^. Moreover, 

r(x) =diag(e(x),Ai(x)) 

$(x) = diag(cr(x),0). 

Here, e and p are symmetric positive definite 3x3 matrices modelling the anisotropic 
permittivity and permeability distribution in the medium, and tr is a symmetric 
positive semi-definite 3x3 matrix which models the anisotropic conductivity 
distribution. In wireless communication, this form can model for example dissipation 
by conductive trees, conductive wires, rainfall, pipes, etc. The operator A can be 
written in block form as 



with 






-83 


82 


^3 





-81 


-82 


81 
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The operators , i — 1,2, 3, can be recovered here from A by putting 

{1 where Am,n = di, 
-1 where A„,„ = -9;, 
elsewhere 

The energy density ^(x, <) is given by 

f(x,t) = l(6(x)|E(x,t)p +Mx)|H(x,t)p). 

The energy flux T{x, t) is described by the Poynting vector 
T{x,t) = E(x,i) X H(x,i). 

5.3. Elastic waves equations 

As already mentioned, also elastic waves can be treated in the general framework of 
symmetric hyperbolic systems. For more details we refer to [331 141j . 

6. The basic spaces and operators 

For our mathematical treatment of time-reversal we introduce the following spaces and 
inner products. We assume that we have J users Uj, j = 1, . . . , J in our system, and in 
addition a base station which consists of K antennas Ak, k = 1, . . . , if. Each user and 
each antenna at the base station can receive, process and emit signals which we denote 
by Sj{t) for a given user Uj, j — 1, . . . , J, and by rk{t) for a given base antenna Ak, 
k = 1, . . . , K . Each of these signals consists of a time-dependent iV- vector indicating 
measured or processed signals of time-length T. In our analysis we will often have to 
consider functions defined on the time interval [T, 2T] instead of [0, T]. For simplicity, 
wc will use the same notation for the function spaces defined on [T, 2T] as we use for 
those defined on [0, T]. It will always be obvious which space we refer to in a given 
situation. 

Lumping together all signals at the users on the one hand, and all signals at the 
base station on the other hand, yields the two fundamental quantities 

s = (si, . . . ,sj) e z 

r = (ri, . . . ,rA') £ Z 

with 

Z = {L2i[0,T]fy, Z = {L2{[0,T]ff. 
The two signal spaces Z and Z introduced above are equipped with the inner products 




The corresponding norms are 

||s||| = (s,s)^, Ilrlll = (r,r)^ 

Each user Uj and each antenna Ak at the base station can send a given signal Sj {t) 
or rfe(i), respectively. This gives rise to a source distribution qj(x, t), j — 1,...,J, 
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or qj.(x, t), k — 1, . . . respectively. Here and in the following we will use in our 
notation the following convention. If one symbol appears in both forms, with and 
without a 'hat' (") on top of this symbol, then all quantities with the 'hat' symbol are 
related to the users, and those without the 'hat' symbol to the antennas at the base 
station. 

Each of the sources created by a user or by a base antenna will appear on the right 
hand side of as a mathematical source function and gives rise to a corresponding 
wave field which satisfies (QJ, lj2l- When solving the system (jSJ, typically certain 
Sobolev spaces need to be employed for the appropriate description of the underlying 
function spaces (see for example lllp. For our purposes, however, it will be sufficient 
to assume that both, source functions and wave fields, are members of the following 
canonical function space U which is defined as 

U = {ue L^in X [0,r])^, u = Oonaf7 X [0,T] , \\u\\u < oo}, 

and which wc have equipped with the usual energy inner product 

(u(t),v(t))y - / / (r(x)u(x,i),v(x,t))^ dxd^, 

J[0,T] Jil 

and the corresponding energy norm 

\Mu = (U> U>(7 • 

Also here, in order to simplify the notation, we will use the same space when 
considering functions in the shifted time interval [T,2T] instead of [0,r]. 

Typically, when a user or an antenna at the base station transforms a signal into a 
source distribution, it is done according to a very specific antenna characteristic which 
takes into account the spatial extension of the user or the antenna. We will model this 
characteristic at the user by the functions 7j(x), j — 1, . . . , J, and for base antennas 
by the functions 7fc(x), k = 1,...,K. With these functions, we can introduce the 
linear 'source operators' Q and Q mapping signals s at the set of users and r at the 
set of base antennas into the corresponding source distributions q(x, t) and q(x, i), 
respectively. They are given as 

J 

Q : Z-^U, q(x,i) ^ Qs^ ^7j(x)s,(i), 

j=i 

K 

Q : Z^U, q(x, t) ^ Qr 7fe(x)rfc(i). 

k=l 

We will assume that the functions 7j(x) are supported on a small neighbourhood 
Vj of the user location d^-, and that the functions 7fc(x) are supported on a small 
neighbourhood Vk of the antenna location dfc. Moreover, all these neighbourhoods are 
strictly disjoint to each other. For example, the functions 7j(x) could be assumed to 
be L2-approximations of the Dirac delta measure S{x — dj) concentrated at the user 
locations d^ , and the functions 7fe(x) could be assumed to be i2-approximations of 
the Dirac delta measure i5(x — dfc) concentrated at the antenna locations dfc. 

Both, users and base antennas can also record incoming fields u £ [/ and 
transform the recorded information into signals. Also here, this is usually done 
according to very specific antenna characteristics of each user and each base antenna. 
For simplicity (and without loss of generality), we will assume that the antenna 
characteristic of a user or base antenna for receiving signals is the same as for 
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transmitting signals, namely jj (x) for the user and 7/j (x) for a base antenna. (The case 
of more general source and measurement operators is discussed in sectional) With 
this, we can define the linear 'measurement operators' M : U ^ Z and M : U ^ Z, 
respectively, which transform incoming fields into measured signals, by 

Sj{t) = {Mu)j = / 7j(x)u(x, t)dx, (j = 1, . . . , J) 
Jn 

rfc(t) = (Mu)fe= /7fe(x)u(x,t)dx, {k^l,...,K) 
Jn 

Finally, we define the linear operator F mapping sources q to states u by 

F : U ->U, = u, 

where u solves the problem (^3), Q- As already mentioned, we assume that the domain 
Vi is chosen sufficiently large and that the boundary dVi is sufficiently far away from 
the users and base antennas, such that there is no energy reaching the boundary in the 
time interval [0, T] (or [T, 2T]) due to the finite speed of signal propagation. Therefore, 
the operator F is well-defined. 

Formally, we can now introduce the two linear communication operators A 
and B which are at the main focus of this paper. They are defined as 

A: Z ^ Z, Av = MFQr, 
B : Z ^ Z, Bs = MFQs. 

The operator A models the following situation. The base station emits the signal r{t) 
which propagates through the complex environment. The users measure the arriving 
wave fields and transform them into measurement signals. The measured signals at 
the set of all users is s{t) — Ar{t). The operator B describes exactly the reversed 
situation. All users emit together the set of signals s{t), which propagate through the 
given complex environment and are received by the base station. The corresponding 
set of measured signals at all antennas of the base station is just r(t) — Bs{t). No 
time-reversal is involved so far. 

7. An inverse problem arising in communication 

In the following, we outline a typical problem arising in communication, which gives 
rise to a mathematically well-defined inverse problem. 

A specified user of the system, say Ui, defines a (typically but not necessarily 
short) pilot signal a{t) which he wants to use as a template for receiving the 
information from the base station. The base station wants to emit a signal r(t) which, 
after having travelled through the complex environment and arriving at the user Ui , 
matches this pilot signal as closely as possible. Neither the base station nor any other 
user except of Ui are required (or expected) to know the correct form of the pilot 
signal a{t) for this problem. As an additional constraint, the base station wants that 
at the other users Z/,-, j > 1, as little energy as possible arrives when communicating 
with the specified user Ui. This is also in the interest of the other users, who want to 
use a different 'channel' for communicating at the same time with the base antenna, 
and want to minimize interference with the communication initiated by user Ui . The 
complex environment itself in which the communication takes place (i.e. the 'channel') 
is assumed to be unknown to all users and to the base station. 
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In order to arrive at a mathematical description of this problem, we define the 
'ideal signal' s{t) received by all users as 

i(0 = («(i),0,...,0)^. (11) 

Each user only knows his own component of this signal, and the base antenna does 
not need to know any component of this ideal signal at all. 

Definition 7.1 The inverse problem of communication: In the terminology 
of inverse problems, the above described scenario defines an inverse source problem, 
which we call for the purpose of this paper the 'inverse problem of communication'. 
The goal is to find a 'source distribution' r{t) at the base station which satisfies the 
'data' s{t) at the users: 

Ar = s. (12) 

The 'state equation' relating sources to data is given by the symmetric hyperbolic 
system (H)), Q. 

Remark 7.1 Notice that the basic operator A in (|12|l is unknown to us since we do 
not know the complicated medium in which the waves propagate. If the operator A 
(together with s) would be known at the base station by some means, the inverse source 
problem formulated above could be solved using classical inverse problems techniques, 
which would be computationally expensive but in principle doable. In our situation, we 
are able to do physical experiments, which amounts to 'applying' the communication 
operator A to a given signal. Determining the operator A explicitly by applying it to 
a set of basis functions of Z would be possible, but again too expensive. We will show 
in the following that, nevertheless, many of the classical solution schemes known from 
inverse problems theory can be applied in our situation even without knowing this 
operator. The basic tool which we will use is a series of time-reversal experiments, 
applied to carefully designed signals at the users and the base station. 

Remark 7.2 A practical template for an iterative scheme for finding an 
optimal signal at the base station can be envisioned as follows. User Ui starts 
the communication process by emitting an initial signal s^i \t) into the complex 
environment. This signal, after having propagated through the complex environment, 
finally arrives at the base station and is received there usually as a relatively long 
and complicated signal due to the multiple scattering events it experienced on its way. 
When the base station receives such a signal, it processes it and sends a new signal 
r^^)(i) back through the environment which is received by all users. After receiving 
this signal, all users become active. The user Ui compares the received signal with 
the pilot signal. If the match is not good enough, it processes the received signal in 
order to optimize the match with the pilot signal when receiving the next iterate from 
the base station. All other users identify the received signal as unwanted noise, and 
process it with the goal to receive in the next iterate from the base station a signal with 
lower amplitude, such that it does not interfere with their own communications. All 
users send now their processed signals, which together define s^^'(i), back to the base 
station. The base station receives them all at the same time, again usually as a long 
and complicated signal, processes this signal and sends a new signal r*^^^ (t) back into 
the environment which ideally will match the desired signals at all users better than 
the previously emitted signal r(^)(t). This iteration stops when all users are satisfied, 
i.e., when user Ui receives a signal which is sufficiently close to the pilot signal a{t), 
and the energy or amplitude of the signals arriving at the other users has decreased 
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enough in order not to disturb their own communications. After this learning process 
of the channel has been completed, the user Ui can now start communicating safely 
with the base antenna using the chosen pilot signal a{t) for decoding the received 
signals. 

Similar schemes have been suggested in 12^1 IS2 EH ^ single-step 

fashion, and in |34[l35j performing multiple steps of the iteration. The main questions 
to be answered are certainly which signals each user and each base antenna needs to 
emit in each step, how these signals need to be processed, at which stage this iteration 
should be terminated, and which optimal solution this scheme is expected to converge 
to. One goal of this paper is to provide a theoretical framework for answering these 
questions by combining basic concepts of inverse problems theory with experimental 
time-reversal techniques. 



8. The basic approach for solving the inverse problem of communication 

A standard approach for solving problem (|12|l in the situation of noisy data is to look 
for the least-squares solution 

ri5 = Minr||Ar-s|i| (13) 

In order to practically find a solution of p3|l . we introduce the cost functional 

J{r) = ^(^Ar - s , Ar - s^^ (14) 

In the basic approach we propose to use the gradient method for finding the 
minimum of H14|l . In this method, in each iteration a correction Sr is sought for a 
guess r which points into the negative gradient direction —A*{Ar — s) of the cost 
functional (|14() . In other words, starting with the initial guess r'-*'-' = 0, the iteration 
of the gradient method goes as follows: 
r(o) = 

r(n-i-i) ^ r(") A* (Ar(") - s), ^^^^ 

where /3'-"-' is the step-size at iteration number n. Notice that the signals r^^"^ are 
measured at the base station, whereas the difference Ar*^"^ — s is determined at the 
users. In particular, s is the pilot signal only known by the user who defined it, 
combined with zero signals at the remaining users, ylr'"^ is the signal received by the 
users at the n-th iteration step. 



9. The adjoint operator A* 

We see from ()15(l that we will have to apply the adjoint operator A* repeatedly 
when implementing the gradient method. In this section we provide practically useful 
expressions for applying this operator to a given element of Z. 

Lemma 9.1 We have 

A* = Q*F*M*. (16) 

Proof: This follows from A ^ MFQ. □ 

Theorem 9.1 We have 

M* = r-ig, M* = r-ig, 

Q* = Air, Q* = MT. ^ ' 
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Proof: The proof is given in appendix A. □ 
Next, we want to find an expression for F*v, v G J7, where F* is the adjoint of 
the operator F. 

Theorem 9.2 Let z be the solution of the adjoint symmetric hyperbolic system 

- - E + - r(x)v(x, t), (18) 

z(x,T) ^ 0, (19) 

z(x,t) = on dnx[0,T]. (20) 

Then 

F*v = T-^{x)z{x,t). (21) 

Proof: The proof is given in appendix B. □ 



Remark 9.1 This procedural characterization of the adjoint operator is often used in 
solution strategies of large scale inverse problems, where it naturally leads to so-called 
'backpropagation strategies'. See for example ^l^^EBEZlEniESI and the references 
given there. 

Remark 9.2 Notice that in the adjoint system H18|) - (|20() 'final value conditions' are 
given at t = r in contrast to Q-Ql where 'initial value conditions' are prescribed at 
t — 0. This corresponds to the fact that time is running backward in (|18|I ~ H20|) and 
forward in O^Q- 



10. The acoustic time-reversal mirror 



In the following we want to define an operator Sa such that F* = T^^SaFSaT holds. 
We will call this operator Sa the 'acoustic time-reversal operator'. We will also 
define the 'acoustic time-reversal mirrors' 7^ and 7^, which act on the signals 
instead of the sources or fields. 

We consider the acoustic system 

+ gradp/(x, i) = qv, (22) 

+ divv/(x,i) =qp, (23) 

v/(x,0) = 0, P/(x,0) ^ in n (24) 

with t £ [0; r] and zero boundary conditions at dfl x [0, T]. We want to calculate the 
action of the adjoint operator F* on a vector {(p, iji)^ E U. 

Theorem 10.1 Let (0,-0)-^ e U and let {va,Pa)'^ be the solution of the adjoint system 
- P-Q^ - gradpQ(x, t) ^ p{x)(f>{x, t) (25) 

-K-^ ~ divva{x,t) = K{x)ilj{x,t) (26) 

v„(x,r) = 0, Pa{x,T) = m (27) 
with t G [O.r] and zero boundary conditions at dQ x [0,T]. Then we have 
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Proof: This theorem is just an appHcation of theorem 19. 21 to the acoustic symmetric 
hyperboUc system. For the convenience of the reader, we will give a direct proof as 
well in appendix C. 

Definition 10.1 We define the acoustic time-reversal operator Sa by putting for 
all q — (qviQp)"'" G U (and similarly for all u = (v,p)"^ ^ U) 

Mc-o^r-*'"!^-") (29) 



qp(x, 2T-t) 

Definition 10.2 We define the acoustic time-reversal mirrors Ta and Ta by 

putting for all r = (rv, rp)^ e Z and all s — (Sv, Sp)^ G Z 

The following lemma is easy to verify. 
Lemma 10.1 We have the following commutations 

MSa=TaM, SaQ ^ QTa, 

MSa=faM, SaQ ^ QTa. ^ ' 

Theorem 10.2 For (0, V)"^ G U we have 

r-5.F5„r(*).f(*) (32) 

Proof: For the proof it is convenient to make the following definition. 

Definition 10.3 Acoustic time-reversal experiment: For given (0, V')^ G U 
define c\y (x, t) ~ p(j) and qp (x, t) = Ktjj and perform the following physical experiment 

d 

p—vtr{x, s) + gradptr(x, s) = - qv(x, 2T ~ s) (33) 

OS 

d 

K—ptr{:i^,s) + diwtr{x,s) =qp(x, 2T-s) (34) 

vtr(x, T) = 0, ptr{yi,T) =0 in Vt, (35) 

with s G [T, 2T] and zero boundary conditions. Doing this experiment means to process 
the data in the following way: Time-reverse all data {(p, ip)^ according to t ^ 2T — s, 
t G [0,T], and, in addition, reverse the directions of the velocities (p — ^ —0. We 
call this experiment the 'acoustic time-reversal experiment'. Notice that the time is 
running forward in this experiment. 

The solution {vtr,Ptr)'^ of this experiment can obviously be represented by 

;::)^ 

In order to show that the so defined time-reversal experiment is correctly modelled by 
the adjoint system derived above, we make the following change in variables: 

T = 2T - S, Vtr = -Vtr, Ptr = Ptr, (37) 

which just corresponds to the application of the operator Sa to {vu- , Ptr)'^ . We have 
T G [0,T]. In these variables the time- reversal system (I33f) - (I35|I gets the form 

d 

^ PQ:^^tr{x, t) - gradptr(x, t) = qv(x, r) (38) 
d 

-K^Ptr(x, t) - divvtr(x, t) =qp(x, t) (39) 
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Vt,(x,T) - 0, Ptr{^,T) =0 in n. (40) 
Taking into account the definition of Qv and qp, we see tliat 

\PtrJ \Pa, 

where Vq and pa solve the adjoint system (|25|l -127 |l . Therefore, according to theorem 

iTTrn 

PtrJ VV" 

Since we have with l|S7|) also 



1 — SaFSa^ 1 / I 1 
PtrJ \V J 

the theorem is proven. □ 



11. The electromagnetic time-reversal mirror 

In the following we want to define an operator Se such that F* ~ T^^SeFSeT holds. 
We will call this operator Se the 'electromagnetic time-reversal operator'. We 
will also define the 'electromagnetic time-reversal mirrors' 7^ and Xj, which act 
on the signals instead of the sources or fields. 
We consider Maxwell's equations 

e-^ - curlH/ + (tE/ = (41) 

M^+curlE/ =qff (42) 

E/(x,0) = 0, H/(x,0) = (43) 

with t e [0, T] and zero boundary conditions. We want to calculate the action of the 
adjoint operator F* on a vector (0,-0)-^ G U . 

Theorem 11.1 Let {(f), ■4>)'^ G U and let (Ea,Ha)-^ he the solution of the adjoint 
system 

(9E 

- e-^ + curlH,(x, t) + aE, = e(x)0(x, t) (44) 

-M-^ - curlEa(x,i) =^(x)?A(x,<) (45) 

E,(x,r) = 0, Ha(x,T)=0 in n (46) 
with < G [0, T] anc? zero boundary conditions. Then we have 

-(:)^(;;:h:U)^ 

Proof: This theorem is just an application of theorem 19.21 to the electromagnetic 
symmetric hyperbolic system. Again, we will give a direct proof in appendix D as 
well. 



Time-reversal and the adjoint method with an application in telecommunication 17 

Definition 11.1 We define the electromagnetic time-reversal operator by 

putting for all q — (q^, q^r)"^ G U (and similarly for all u = (E, H)'^ ^ U) 

(5eq)(x,<) = ( , _ (48) 



qff(x, 2r-t) 

Definition 11.2 We define the electromagnetic time-reversal mirrors %. and 

Te by putting for all r = (r^;, th)^ & Z and for all s — {se, sh)^ 6 Z 

The following lemma is easy to verify. 
Lemma 11.1 We have the following commutations 

MSe = TeM, SeQ = QTe, 

MSe ^ feM, SeQ^Qfe. ^ ' 

Theorem 11.2 For {(j),ip)^ G U we have 

r-.5.F..r(^t)^^-(:) <-) 

Proof: For the proof it is convenient to make the following definition. 

Definition 11.3 Electromagnetic time-reversal experiment: For a given vector 
{(pjtp)^ e U define q£;(x, i) = ecj) and q/f(x, i) = piip and perform the physical 
experiment 

e^Et,(x,s) - curlHt,(x,s) + (7Et,(x,s) = -qB(x,2r-s) (52) 

OS 

^|-Hi,(x, s) + curlEt,(x, s) = qni^, 2T - s) (53) 

OS 

Etr(x,r) = 0, Ht^(x,r) = in Q, (54) 

with s e [T, 2T] and zero boundary conditions. Doing this experiment means to process 
the data in the following way: Time-reverse all data according to t ^ 2T — s, t G [0,T], 
and, in addition, reverse the directions of the electric field component by (j) ^ —(p. We 
call this experiment the 'electromagnetic time-reversal experiment'. Notice that the 
time is running forward in this experiment. 

The solution (EtrjUtr)'^ of this experiment can obviously be represented by 



In order to show that the so defined time-reversal experiment is correctly modelled by 
the adjoint system derived above, we make the following change in variables: 

r = 2T - s, Etr - -Etr, B.tr = ntr, (56) 

which just corresponds to the application of the operator Se to (Ef^, Ht^)^- We have 
r G [0,T]. In these variables the time-reversal system (I52f) - (|54|l gets the form 

d . 

- e— Etr(x,T) + curlHtr(x, r) + aEtr{x,s) = qE(x, r), (57) 

OT 

d - 

-//— Htr(x, r) - curlEtr(x,T) =qH(x, r) (58) 
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= TF 

Since we have with H55|l . H5fci|) also 

= SeFSeT 




Et^(x,r) = 0, Ht^(x,r) in Q. (59) 
Taking into account the definition of and qn, we see that 
/ Et, \ Ea 

where Eg and He solve the adjoint system (|44|I - H46|I . Therefore, according to theorem 

irm 



the theorem is proven. □ 

Remark 11.1 For electromagnetic waves, there is formally an alternative way 
to define the electromagnetic time-reversal operator, namely putting for all q = 

qB(x, 2T-t) 
-qH(x, 2T-t), 

accompanied by the analogous definitions for the electromagnetic time-reversal 
mirrors. With these alternative definitions, theorem 111.21 holds true as well, with only 
very few changes in the proof. Which form to use depends mainly on the preferred 
form for modelling applied antenna signals in the given antenna system. The first 
formulation directly works with applied electric currents, whereas the second form is 
useful for example for magnetic dipole sources. 

12. Time-reversal and the adjoint operator A* 

Define S — Sa, T — Ta, T — % ioi the acoustic case, and S = Se, T = T — % 
for the electromagnetic case. We call S the time-reversal operator and T, T the 
time-reversal mirrors. We combine the results of lemma IIU. II and lemma [Tl.ll into 
the following lemma. 

Lemma 12.1 We have the following commutations 
MS = TM, SQ = QT, 

MS = TM, SQ = QT. ^ ' 

Moreover, combining theorem 1 1 (J . 21 and theorem 1 11. 21 we get 
Theorem 12.1 For {(j),ip)^ G U we have 

V-^SFST i'^ ] ^ F* i'^ ] . (61) 



(5eq)(x,<) = 



With this, we can prove the following theorem which provides the fundamental link 
between time-reversal and inverse problems. 
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Theorem 12.2 We have 

A* = TBT. (62) 

Proof: Recall that the adjoint operator A* can be decomposed as A* = Q*F*M*. 
With theorem l9.1l theorem 1 12. II and lemma [T2. II it follows therefore that 

A* = Mrr~^SFsrr-^Q 

= MSFSQ 
= TMFQT 
= TBT, 

which proves the theorem. □ 

Remark 12.1 The above theorem provides a direct link between the adjoint operator 
A*, which plays a central role in the theory of inverse problems, and a physical 
experiment modelled by B. The expression TBT defines a 'time-reversal experiment'. 
We will demonstrate in the following sections how wc can make use of this relationship 
in order to solve the inverse problem of communication by a series of physical time- 
reversal experiments. 

Remark 12.2 We only mention here that the above results hold as well for elastic 
waves with a suitable definition of the elastic time-reversal mirrors. We leave 
out the details for brevity. 



13. Iterative time-reversal for the gradient method 

13.1. The basic version 

The results achieved above give rise to the following experimental procedure for 
applying the gradient method (|15(l to the inverse problem of communication as 
formulated in sections and |H1 First, the pilot signal s{t) is defined by user Ui 
as described in (|ll|l . Moreover, we assume that the first guess r^^\t) at the base 
station is chosen to be zero. Then, using theorem 112.21 we can write the gradient 
method H15|l in the equivalent form 

r(0) = 

j.(n+i) = r(") + /?(")rsr(s- Ar(")), ^""^^ 

or, expanding it, 

r(0) = 

s(") = s-Ar(") (64) 

j.(n+l) _ j.(ri) _j_ ^(")7'S'7's(") 

In a more detailed form, we arrive at the following experimental procedure for 
implementing the gradient method, where we fix in this description /j'"^ = 1 for 
all iteration numbers n for simplicity. 

(i) The user Ui chooses a pilot signal a{t) which he wants to use for communicating 
with the base station. The objective signal at all users is then s{t) = 
(a(t),0,...,0)^. The initial guess Y^^>{t) at the base station is defined to be 
zero, such that s^*'-' = s. 

(ii) The user lAi initiates the communication by sending the time-reversed pilot signal 
into the environment. This signal is Ts'^'^\t). All other users are quiet. 
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(iii) The base station receives the pilot signal as BTs'^^^ (t). It time-reverses this signal 
and sends this time-reversed form, namely r(^)(t) ~ TBTs'^^^t), back into the 
medium. 

(iv) The new signal arrives at all users as Ar^-^'>{t). All users compare the received 
signals with their components of the objective signal s{t). They take the difference 
s^^\t) = s(t) — Ar'^^^t), and time-reverse it. They send this new signal Ts(^)(t) 
back into the medium. 

(v) The base station receives this new signal, time-reverses it, adds it to the previous 
signal r(^)(t), and sends the sum back into the medium as r^^\t). 

(vi) This iteration is continued until all users are satisfied with the match between the 
received signal Ar'-^-^t) and the objective signal s{t) at some iteration number 
n. Alternatively, a fixed iteration number n can be specified a-priori for stopping 
the iteration. 

Needless to say that, in practical implementations, the laboratory time needs to be 
reset to zero after each time-reversal step. 

Remark 13.1 The experimental procedure which is described above is practically 
equivalent to the experimental procedure which was suggested and experimentally 
verified in j34ll35| . Therefore, our basic scheme provides an alternative derivation and 
interpretation of this experimental procedure. 

Remark 13.2 We mention that several refinements of this scheme are possible and 
straightforward. For example, a weighted inner product can be introduced for the 
user signal space Z which puts different preferences on the satisfaction of the user 
objectives during the iterative optimization process. For example, if the 'importance' 
of suppressing interferences with other users is valued higher than to get an optimal 
signal quality at the specified user Ui , a higher weight can be put into the inner product 
at those users which did not start the communication process. A user who does not 
care about these interferences, simply puts a very small weight into his component of 
the inner product of Z. 

Remark 13.3 Notice that there is no mechanism directly built into this procedure 
which prevents the energy emitted by the base antenna to increase more than the 
communication system can support. For example, if the subspace of signals 

Zo := {r{t) : Ar = 0} 

is not empty, then it might happen that during the iteration described above (e.g. 
due to noise) an increasing amount of energy is put into signals emitted by the base 
station which are in this subspace and which all produce zero contributions to the 
measurements at all users. More generally, elements of the subspace of signals 

Z, {r{t) : \\Ar\\^ < e\\r\\z}, 

for a very small threshold < e << 1, might cause problems during the iteration if 
the pilot signal s{t) chosen by the user has contributions in the subspace AZ^ (i.e. in 
the space of all s = Ar with r G Z^). This is so because in the effort of decreasing the 
mismatch between Ar and s{t), the base antenna might need to put signals with high 
energy into the system in order to get only small improvements in the signal match 
at the user side. Since the environment (and therefore the operator A) is unknown 
a-priori, it is difficult to avoid the existence of such contributions in the pilot signal. 



Time-reversal and the adjoint method with an application in telecommunication 21 

One possible way to prevent the energy emitted by the base station to increase 
artificially would be to project the signals r^") (t) onto the orthogonal complements of 
the subspaces Zq or (if they are known or can be constructed by some means) prior 
to their emission. Alternatively, the iteration can be stopped at an early stage before 
these unwanted contributions start to build up. (This in fact has been suggested in 

ISIESl) 

In the following subsection we introduce an alternative way of ensuring that the 
energy emitted by the base station stays reasonably bounded in the effort of fitting 
the pilot signal at the users. 



13.2. The regularized version 

Consider the regularized problem 

TLSr = Minr(||Ar - s|l| + A|lr|||) (65) 

with some suitably chosen regularization parameter A > 0. In this problem formulation 
a trade-off is sought between a signal fit at the user side and a minimized energy 
emission at the base station. The trade-off parameter is the regularization parameter 
A. Instead of (|14|l we need to consider now 

J(r) i{|^r-s, Ar-sJ)^ + ^{Jr, r^^ (66) 

The negative gradient direction is now given by —A*{Ar — s) — Ar, such that the 
regularized iteration reads: 

r(0) ^ 



j.(n+l) ^ 



[.(") + /3(") (TBf{s - Ar(")) - Ar(")) , 



(67) 



where we have replaced A* by TBT . The time-reversal iteration can be expanded 
into the following practical scheme 

r(0) = 

s(") = s-Ar(") (68) 
j.(n+i) _ j.(n) _|_ ^(")7'S7's(") — /3(")Ar("). 

Comparing with (|64|) . we see that the adaptations which need to be applied in the 
practical implementation for stabilizing the basic algorithm can easily be done. 



14. The minimum norm solution approach 

14-. 1. The basic version 

In this section we want to propose an alternative scheme for solving the inverse problem 
of communication. As mentioned above, a major drawback of the basic approach (|13|l 
is that the energy emitted by the base station is not limited explicitly when solving 
the optimization problem. The regularized version presented above alleviates this 
problem. However, we want to mention here that, under certain assumptions, there 
is an alternative scheme which can be employed instead and which has an energy 
constraint directly built in. Under the formal assumption that there exists at least 
one (and presumably more than one) solution of the inverse problem at hand (i.e. the 
'formally underdetermined case'), we can look for the minimum norm solution 

Miurllrll^ subject to Ar = s. (69) 
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In Hilbert spaces this solution has an exphcit form. It is 

VMN = A*{AA*)-^~s. (70) 

Here, the operator (^4^4*)"^ acts as a filter on the pilot signal s. Instead of sending 
the pilot signal to the base station, the users send the filtered version of it. Certainly, 
a method must be found in order to apply the filter {AA*)~^ to the pilot signal. One 
possibility of doing so would be to try to determine the operator AA* explicitly by a 
series of time- reversal experiments on some set of basis functions of Z , and then invert 
this operator numerically. However, this might not be practical in many situations. (It 
certainly would be slow and it would involve a significant amount of signal-processing, 
which we want to avoid here.) Therefore, we propose an alternative procedure. First, 
we notice that there is no need to determine the whole operator {AA*)~^ , but that 
we only have to apply it to one specific signal, namely s. Let us introduce the short 
notation 

C = AA*. 

In this notation, we are looking for a signal s € Z such that Cs = s. We propose to 
solve this equation in the least squares sense: 

s = Mins||Cs-s|||. (71) 

Moreover, as a suitable method for practically finding this solution, we want to use 
the gradient method. Starting with the initial guess sq = 0, the gradient method 
reads 

s(°) = 

s(n+i) = s(") - /3(")C*(Cs(") -s), ^^^^ 

where C* is the adjoint operator to C and /3^"^ is again some step-size. Expanding 
this expression, and taking into account C* = C and A* = TBT, we arrive at 
s(o) = 

s(n-M) ^ s(") + /3("'ATBr(s-ylT.Bfs(")). ^'^^^ 

In the practical implementation, we arrive at the following iterative scheme: 

Initialize gradient iteration: s'^^ = 
for n = 0, 1, 2, . . . do: 

s("+5) = s-Ar("+5) 

g(ri+l) _ g(n) ^ p(n) j^j.(n+l) 

Terminate iteration at step n: s = s^"^ 

Final result: tmjv = TBTs. 

This iteration can be implemented by a series of time-reversal experiments, without 
the need of heavy signal-processing. The final step of the above algorithm amounts 
to applying A* to the result of the gradient iteration for calculating s — {AA*)^^s, 
which yields then tmn- This will then be the signal to be applied by the base station 
during the communication process with the user Ui . 



Time-reversal and the adjoint method with an application in telecommunication 23 



14-.2. The regularized version 

In some situations it might be expected that the operator C is iU-conditioncd, such that 
its inversion might cause instabihties, in particular when noisy signals are involved. 
For those situations, a regularized form of the minimum norm solution is available, 
namely 

rMNr = A*{AA* + Xl^y^s (75) 

where denotes the identity operator in Z and A > is some suitably chosen 
regularization parameter. The necessary adjustments in the gradient iteration for 
applying {AA* + XI ^)^^ to s are easily done. We only mention here the resulting 
procedure for the implementation of this gradient method by a series of time-reversal 
experiments: 

Initialize gradient iteration: s^°) = 
for n = 0, 1, 2, . . . do: 

j.(n+i) ^ TSTs("' 

j.(n+l) ^ TSts("+5) 

g(n+i) _ s(") + /3(")Ar("+^) + /3("^s("+5). 
Terminate iteration at step n: s ~ s*^"^ 

Final result: ymnt — TBTs. 

Again, the last step shown above is a final application of A* to the result of the 
gradient iteration for calculating s = [AA* + A/^)^^s, which yields then ymnt- This 
will then be the signal to be applied by the base station during the communication 
process with the user Ui . 



15. The regularized least squares solution revisited 

We have introduced above the regularized least squares solution of the inverse problem 
of communication, namely 

TLSr = Minr(||Ar - + A||r|||) (77) 

with A > being the regularization parameter. In Hilbert spaces, the solution of (|77|l 
has an explicit form. It is 

TLSr ^ iA*A + XIzr'A*s, (78) 

where Iz is the identity operator in Z. It is therefore tempting to try to implement also 
this direct form as a series of time-reversal experiments and compare its performance 
with the gradient method as it was described above. As our last strategy which 
we present in this paper, we want to show here that such an alternative direct 
implementation of (|77|) is in fact possible. 

Notice that in H77|) the filtering operator {A* A + XI z)^^ is applied at the base 
station, in contrast to the previous case where the user signal was filtered by the 
operator {AA* + Xl2)~^ ■ Analogously to the previous case, we need to find a practical 
way to apply this filter to a signal at the base station. We propose again to solve the 
equation 

iA*A + XIz)r = i (79) 
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(80) 



in the least squares sense, where f = A*s. Defining 
C = A*A + \Iz, 

and using C* — C and A* ~ TBT, we arrive at the following gradient iteration for 
solving problem (|79|l : 

r(0) = 

i.(n+i) ^ r(") + /?(")(rsf A + A/z)(f - (rBrA + A/z)r(") 

This gives rise to the following practical implementation by a series of time-reversal 
experiments: 

User sends pilot signal to base station: f = TBTs 

Initialize gradient iteration: r^"^ = 
for n = 0, 1, 2, . . . do: 

g(n+5) _ Ar^"^ 

r(«+3) = f _ rBrs("+^) - Ar(") (81) 
g(n+i) ^ Ar^"+i^ 

j.(n+i) _ r(") + /3(")Ti?Ts("+^) + /3(")Ar("+5). 
Terminate iteration at step h: r — r'"^ 
Final result: J^LSr = f . 
TLSr will then be the signal to be applied by the base station during the communication 
process with the user Ui . 

16. Partial and generalized measurements 

In many practical applications, only partial measurements of the whole wave-field are 
available. For example, in ocean acoustics often only pressure is measured, whereas 
the velocity field is not part of the measurement process. Similarly, in wireless 
communication only one or two components of the electric field might be measured 
simultaneously, but the remaining electric components and all magnetic components 
are missing. We want to demonstrate in this section that all results presented above 
are valid also in this situation of partial measurements, with the suitable adaptations. 

Mathematically, the measurement operator needs to be adapted for the situation 
of partial measurements. Let us concentrate here on the special situation that only one 
component u^, (^ S {1, 2, 3, ...}) of the incoming wave field u is measured by the users 
and the base station. All other possible situations will then just be combinations of 
this particular case. It might also occur the situation that users can measure a different 
partial set of components than the base station. That case also follows directly from 
this canonical situation. 

We introduce the new signal space at the base station Y = {L2[0,T])^ and the 
corresponding 'signal projection operator' P^, by putting 

: Z^Y , P,{ri{t),...,rK{t)f = v,{t). 

We see immediately that its adjoint P* is given by 

P::Y^Z , P;r,(t) = (0,...,0,r,(i),0,...,0)^ 

where v^{t) appears on the right hand side at the i^-th position. Our new measurement 
operator M^, and the new source operator Q^, are then defined by 

M^:U^Y, AUu = P^Mu 

Q,:Y^U, Q,v, = QP:t,. ^^^> 
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Analogous definitions are done for Y, P^,, My and Qy at the users. 

Obviously, we will have to replace now in the above derivation of the iterative 
time- reversal procedure all measurement operators M by My (and M by My) and 
all source operators Q by Qy (and Q by Qu)- In particular, the new 'communication 
operators' are now given by 

A.y '. Y > Y. AyVy — MyFQyVy, 

By-.Y^Y, ByS^ = MyFQySy. ^ ' 

In the following two theorems we show that the main results of this paper carry over 
to these newly defined operators. 

Theorem 16.1 We have 

m; = r-ig„ m; = r-iQ„ 

Ql = MyT, Ql = MyT. ^ ' 

Proof: The proof is an easy exercise using (|82|l and theorem 19.11 □ 

Theorem 16.2 It is 

A*y = TByf. (85) 

Proof: The proof is now identical to the proof of theorem ll2.2l using theorem ll6.ll 
instead of theorem 19. II □ 

Remark 16.1 In fact, it is easy to verify that all results of this paper remain valid 
for arbitrarily defined linear measurement operators 

Mu -.U^Zu, Ma-.U^ Za, 

where Zu and Za are any meaningful signal spaces at the users and the base antennas, 
respectively. The only requirement is that it is experimentally possible to apply signals 
according to the source operators defined by 

Qu ■■ Zu^ U, Qu - TMli 

Qa- Za^ U, Qa = TMX 

where Mj^ and M^ are the formal adjoint operators to Mu and Ma with respect to the 
chosen signal spaces Zu and Za- In addition, the measurement and source operators 
as defined above are required to satisfy the commutation relations as stated in lemma 
112.11 Under these assumptions, we define the generalized communication operators A 
and B by 

A: Za^Zu, Ar = MuFQat, 
B: Zu^ Za, Bs = MaFQus. 

Now the proof to theorem 1 1 2 . 21 directlv carries over to this generalized situation, such 
that we have also here 

A* ^ TBT. 

This yields iterative time-reversal schemes completely analogous to those presented 
above. 
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17. Summary and future directions 

We have derived in tliis paper a direct link between the time-reversal technique and 
the adjoint method for imaging. Using this relationship, we have constructed several 
iterative time-reversal schemes for solving an inverse problem which arises in ocean 
acoustic and wireless communication. Each of these schemes can be realized physically 
as a series of time-reversal experiments, without the use of heavy signal processing or 
computations. One of the schemes which we have derived (and which we call the 'basic 
scheme'), is practically equivalent to a technique introduced earlier in |34[ I35| using 
different tools. Therefore, we have given an alternative theoretical derivation of that 
technique, with a different mathematical interpretation. The other schemes which we 
have introduced are new in this application. They represent either generalizations of 
the basic scheme, or alternatives which follow different objectives. 

Many questions related to these and similar iterative time-reversal approaches 
for telecommunication are still open. The experimental implementation has been 
investigated so far only for one of these techniques in |34[ I35j . for the situation 
of underwater sound propagation. A thorough experimental (or, alternatively, 
numerical) verification of the other schemes is necessary for their practical evaluation. 
An interesting and practically important problem is the derivation of quantitative 
estimates for the expected focusing quality of each of these schemes, for example 
following the ideas of the work performed for a single step in [S]. Certainly, it is 
expected that these estimates will again strongly depend on the randomness of the 
medium, on the geometry and distribution of users and base antennas, and on technical 
constraints as for example partial measurements. Also, different types of noise in the 
communication system need to be taken into account. The performance in a time- 
varying environment is another interesting issue of practical importance. All schemes 
presented here can be adapted in principle to a dynamic environment by re-adjusting 
the constructed optimal signals periodically. Practical ways of doing so need to be 
explored. 

Finally, we want to mention that speed is an important factor in the application of 
communication. In the time-reversal experiment, in a certain sense a physical process 
is used for applying a mathematical operator (namely A*) in order to solve a given 
inverse problem. In a complex environment, this can be much faster than using a 
powerful computer, even if the complex environment would be known to the users and 
the base station. 
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Appendix A: Proof of theorem 19.11 
For arbitrary r d Z and u G C/ we liave 

(Mu,r}^ - / / / 7fc(x)u(x,i)dx, r{t)\ dt 

^.^^ J[Q,T\ \Jn I N 

= [ [ (r(x)u(x,i), r-i(x)^7fc(x)r(i)\ d^dt. 
J\o.T] Jn \ , / .. 



Therefore, 



with 



K 



fe=l / AT 



M*r = r-i(x)^7fe(x)r(i) = r-i(x)Qr. 
fc=i 

Doing the analogous calculation for M* we get 

M* = r-^Q, M* = r-^Q. (A.l) 

Taking the adjoint of IjA.ip we see that 

Q* ^ MT, Q* = Mr 
holds as well. This proves the theorem. □ 

Appendix B: Proof of theorem 19.21 

We have the following version of Green 's formula 

+ / / (r(x)u(x,t),v(x,i))^ dxrft + / / (q(x,t), z(x,i))jY dxdt 
J[O.T]Jfl J[o,T]Jn 



+ / / (r(x)u(x,t),v(x,t))^ dxdi + / / (q(x,t), z(x,i))jY dxdt 
J[(),T]Jn J[o,T]Jn 

+ f (r(x)u(x,r),z(x,T))j^ dx - / (r(x)u(x,0),z(x,0)}^ dx 
Jn Jn 

+ / / V(i^'u,z)^i/,(x)dadt, 

J[a,T] Jon ^^-^ 

where n(x) — (i'i(x), j^2(x), 1^3 (x)) is the outward normal at dO. in the point x. Notice 
that we have augmented Green's formula in ljA.2p by some terms which appear in 
identical form on the left hand side and on the right hand side. 
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We will assume here that the boundary is far away from the sources and receivers 
and that no energy enters fl from the outside, such that during the time interval of 
interest [0, T] all fields along this boundary are identically zero. This is expressed by 
the boundary conditions given in Q and jSUJ. Let u(x, t) be a solution of (Q),©)®, 
and z(x, t) a solution of (|18() - (|20|l . Then the first term on the left hand side of (|A.2p 
and the third term on the right hand side cancel each other because of . The second 
term on the left hand side and the first term on the right hand side cancel each other 
because of H18|l . The {t — r)-term and the (t = 0) term vanish due to l|19|l and 

respectively, and the boundary integral vanishes because of Q and H2UI) . The 
remaining terms (i.e. the third term on the left hand side and the second term on the 
right hand side) can be written as 

with F*v = r-i(x)z(x, t) as defined in □ 
Appendix C: Direct proof of theorem 110.11 



We prove the lemma by using Greens formula: 



T 

Jfl 



dvf dpf 



dyidt (A.3) 



/ / [pvf(f) + Kpf^] d:x.dt+ / / [qvVa + qpPa] dxdt 
Jo Jn Jo Jn 



Jn 

T 



d-Va ,. dpa , 

'^fP~Ql P/QivVq - Pfi^-g^ v/gradpa 



dx.dt 



+ / [pvfcj) + npfil^] dxdt+ / / [qvVa + qpPa] dxdt 
Jo Jn Jo Jn 

+ / {va-n)pfdadt + / / f ■ n)pa dadt 
Jo Jon Jo Jon 

+ I P[(v/Va)(x,T)-(v/V,)(x,0)]dx + / K.[{pfPa){x,T) - {pfPa){x,Q)\dx. 

Jn Jn 
This equation has the form ljA.2|) . Notice that we have augmented Green's formula 
in ljA.3|) . as already shown in ljA.2|) . by some terms which appear in identical form on 
the left hand side and on the right hand side. 

The first term on the left hand side of equation ljA.3|l and the third term on the 
right hand side cancel each other due to H22I) . H23|1 . The second term on the left hand 
side and the first term on the right hand side cancel each other because of (|25|) . (|26|) . 
The {t — T)-terms and the {t — 0)-terms vanish due to (|27() . H24|l . respectively, and 
the boundary terms vanish because of zero boundary conditions. We are left over with 
the equation 

/ / + npfi)] dxdt = / [qvVa + qpPa] dxdt 

Jo Jn Jo Jn 

Defining F* by (|28|) . this can be written as 



Therefore, F* is in fact the adjoint of F, and the lemma is proven. □ 
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Appendix D: Direct proof of theorem lll.il 



We prove the lemma by using Green 's formula: 



Jn 



e^Ea - curlH/Ea + aEfEa 



dt 



dt 



Ha + curlE/Ha 



d'x.dt 



[ [ [eE/0 + AiH/V^] dxdt + [ [ [ciEEa + q^Ha] dxdt 
Jo Jn Jo Jn 



Ete 



Jn I 

T 



dEg 

dt 



H/curlEa + EfoEa - Hj^J. 



dUg 

dt 



E/curlHa 



dxdt 



[ [ [eE/0 + /xH/V'] dxdt + [ [ [qsEg + q^Ha] dxdi 
Jo Jn Jo Jn 

T 



T 

Jn 



Jan 



Eg xHf ■ n dadt 



Jan 



Ef X Hg ■ n dadt 



(A.4) 



+ 



/ e[{EfEg){T)-{EfEg){0)]dx+ f A^[(HyHa)(r) - (H^Ha)(0)] dx 
Jn Jn 



This equation has the form ljA.21) . Notice that we have augmented Green's formula 
in ljA.4|) . as already shown in IIA.21) . by some terms which appear in identical form on 
the left hand side and on the right hand side. 

The first term on the left hand side of equation HA.4|I and the third term on the 
right hand side cancel each other because of H41|) and 142|) . The second term on the 
left hand side and the first term on the right hand side cancel each other because of 
(O, (gSJi. The {t = 0)-terms and the {t = r)-terms vanish due to l|33Il and gSl- The 
boundary terms vanish because of zero boundary conditions. We are left over with 
the equation 

/ / [eEf(f> + nHfip] dTidt = / / [q_EEa + q^Ha] d:x.dt. 
Jo Jn Jo Jn 

Defining F* by (|47|) . this can be written as 



F 



f qs 

\^.H 



f qfi 
vqH 



, F* 



Therefore, F* is in fact the adjoint of F , and the lemma is proven. 



□ 
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